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1. Introduction 


The estimation of non-perforation probabilities is of significant interest to the developers of 
annor systems, whose objective is to provide protection for vehicles and personnel against 
enemy threats. In particular, the following questions may arise: 

1. What is the probability that homogenous annor of a given thickness will not be 
perforated by a specific enemy threat? 

2. How thick should the armor be if one is to be highly confident that the probability of 
non-perforation with this threat is high? 

Conversely, these probabilities are of interest to the developers of armor-piercing projectiles, 
whose objective is to perforate enemy annor. From this point of view, one might ask 

3. What is the probability that a projectile will be able to perforate an enemy armor of a 
given thickness? 

For a particular projectile, non-perforation 1 occurs if its penetration depth, denoted by x, is less 
than the thickness of the armor, denoted by x 0 . Now consider the random variable X, which 

represents the penetration depth of a randomly selected projectile. The probability of non¬ 
perforation is an unknown, fixed constant within the interval [0, 1] and is written as P(X < x 0 ). 

An estimate for the true non-perforation probability will be based on sample data and as such, is 
subject to sampling variation. However, a point estimate for the non-perforation probability is a 
single value and therefore does not give any sense of the uncertainty associated with it. 

Preferable to a point estimate is an interval estimate that contains the non-perforation probability 
with high probability, i.e., a confidence interval. The width of the interval provides a 
quantitative measure of the estimation error. 

The problem of estimating P(X < x n ) has its origins in acceptance sampling plans for statistical 

quality control in the 1950’s. Liebennan and Resnikoff (1955) and Barton (1961) yielded 
expressions for the minimum variance estimator. It was not until the 1970’s that confidence 
interval estimators surfaced. In their 1977 paper, Owen and Hua derived lower confidence limits 
for P{X > x 0 ) and described how to obtain one-sided and two-sided limits for P(X <x 0 ). 

However, their derivation contains several errors in notation and as a result, is quite difficult to 
follow. Owen and Hua’s paper also provided tables that assist in the calculation of lower 
confidence limits for P(x > x 0 ); these tables were dramatically expanded in Odeh and Owen 


1 For the remainder of this report, we only discuss non-perforation events, their probabilities of occurrence, and 
interval estimates for these probabilities. Analogous results for the probability of perforation, P(X > x 0 ), are easily 
derived if we use the axiomatic property P(X > x 0 ) = l — P(X < x 0 ). 
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(1980, tables 7.1 through 7.7). Odeh and Owen’s tables are cited in the more recent engineering 
statistics text by Hahn and Meeker (1991, section 4.5). 

In this report, the recent concept of generalized inference is used to develop one- and two-sided 
confidence intervals for the probability of non-perforation when it can be assumed that the 
penetrations follow a normal distribution. Although this may be seen as an attempt to “reinvent 
the wheel,” it is meant to show the wide applicability of the generalized approach to confidence 
interval construction. As pointed out by Weerahandi (2004, section 1.7), generalized inference 
may be used to obtain confidence intervals for any function of the normal distribution 
parameters. 

Solutions for the limits of the intervals require iterative calculation of percentiles of non-central t 
distributions. As an alternative to the tables by Odeh and Owen and to avoid the interpolation 
errors inherent in their usage, software code is provided in the appendices so that the practitioner 
can easily obtain exact solutions for the interval limits. 


2. General Solution 


Suppose that n depths of penetration x v x 2 ,...,x n are collected in an experiment for which the 

observations are assumed to come from a nonnal distribution with mean p and standard deviation 
a. Historically, this assumption has been found to be tenable when the projectiles are of the same 
ammunition type with (approximately) equal muzzle velocities, obliquities, etc., and the target is a 
homogeneous armor. For the purpose of this report, we assume that the armor is thick enough to 
guarantee that no projectiles will perforate. From this sample, an estimate is desired for the 
probability of non-perforation, i.e., the probability that a future penetration, X, from the same 
population of projectiles will be less than or equal to some specified depth, x 0 . This probability is 

written as P(X < x 0 ) and is known as the cumulative distribution function (CDF) evaluated at x 0 , 
or F x (x 0 ). Notice that 


P(X<x 0 ) = P 


X ~M < *o ~M 


= P\ 


Z< 


x 0 -/j 


= <J> 


) 


A) 


V 


) 


( 1 ) 


in which Z is a standard normal random variable and d>() is the standard normal cumulative 

distribution function (CDF) commonly tabulated in statistics texts and available in many 
software packages. Therefore, calculating a simple point estimate for P(X < x () ) is 


straightforward: one computes the standard normal CDF at an estimate of 


x 0 -fi 


. We obtain 


1 " 

this estimate by substituting the sample mean (i = -Vi j ) and the sample standard deviation 

n ;=1 


2 



(s 


-XU, ~ A') 2 ) for the population mean ( fi ) and the population standard deviation (<x), 

n~ 1 m 


respectively 2 . That is, 


F x (x a ) = P(X<x„)=P 


Z< X «ZA 




X n -x 
Z<— - 


$ 


^ x 0 — x^ 


V a ) k s ) \ s j 

However, any simple point estimate for F x (jc 0 ) does not give information about its precision. 

What we would prefer to construct is a confidence interval or perhaps a confidence bound for 
F x (x 0 ) that by its width gives some sense of the degree of uncertainty associated with it. 3 


( 2 ) 


We begin by recalling from equation 1 that F x (x 0 ) = <T> 


x 0 // ^ 


V o ) 

monotonic function, we can easily attain a confidence interval for <f> 


. Since <f>( ) is a strictly 
by deriving a 


f x 0 ~M^ 


confidence interval for 


x 0 ~fi 


k o- J 

and then applying the function <J>( ) to the resulting confidence 


limits (see Mood, Graybill, and Boes, 1974, page 378). 


X — LI 

The following derivation of a confidence interval for 0 = ——— is based on a method by 

(7 

Weerahandi (1995) in which he uses a generalized pivotal quantity to obtain a confidence 

interval for the coefficient of variation, —, for a normal population with mean p and standard 

M 

deviation a. The construction of a generalized pivotal quantity requires one to find a function R 
with arguments 

1. X and S 2 (independent and sufficient statistics for the random sample X l ,..., X n ), 

2. x and 5 2 (the observed values of X and S 2 ), 

3. 6 , the parameter of interest, and 

4. perhaps a vector, £, of additional unknown (nuisance) parameters. 


2 Although this estimate is easily obtained, it does not have the desirable property of being unbiased. An unbiased 
estimator with minimum variance is presented in either Lieberman and Resnikoff (1955) or Barton (1961). 


3 Note the distinction between the problems of solving for a confidence interval for <t> 


Xq~F 


and solving for a 


tolerance interval. With the former, the percentile value ( X 0 ) is specified and an interval for its associated CDF 

*o~M ' 


value, <f> 


, is desired. On the other hand, with tolerance interval construction, the CDF value (p) is 


specified while an interval for the associated percentile, p + z- <7 , is desired. 


3 



This function, denoted in full by r{x, S 2 ; x, s 2 , 0 , g), must satisfy the following two conditions: 

Condition 1: r[x, S 2 ; x, s 2 ,6, £) has a distribution that is free of any unknown parameters. 

Condition 2: The observed value of R, i.e., r = r{x, s 2 ; x, s 2 , 0, £), is equal to 6. 

If such a function can be found, then it is a generalized pivotal quantity for 6 , and its percentiles 
can be used to obtain confidence intervals for 0 . Since their introduction in the 1980’s, 
generalized pivotal quantities have been regarded as quite challenging to derive. Even Weerahandi 
(1993) states “... the construction of pivotals requires some intuition.” In essence, one obtains the 
pivotal quantity by working backwards from the expression for 0 . The mathematical “tricks” of 
adding 0 and/or multiplying by 1 are adroitly employed to link all unknown parameters in 6 with 
functions of the sufficient statistics whose distributions are parameter free. Once the linking is 
complete, all remaining random variables are converted to their observed values. 

Step 1 : The sufficient statistics for a normally distributed random sample are X and S 2 . 

Random variables based on these statistics whose distributions are free of unknown parameters 
include 

X — LI 1 

1. Y =-, distributed as a normal random variable with mean zero and variance — ; and 

<7 n 

S 2 

2. V = (n -1) —-, distributed as a chi-square random variable with n - 1 degrees of freedom. 

<J~ 


Step 2: We attempt to construct a random variable from 6 that involves Y and/or V. First, 
notice that by adding zero in the form of [-X + x) to the numerator, 0 can be rewritten in the 
following manner: 

e _ x 0 -ju _ x 0 -X + X->u _ x 0 -X | X-ju 

a G G G 


X — LI 

Now the parameter ji is linked to the random variable-which was previously noted as 

G 

X — u 

having a parameter-free distribution. When we make the substitution Y =-—, the parameter 


// is removed from 6 : 


0 = 


-X 


+ Y. 


G 


(4) 
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It remains to manipulate the left addend of equation 4 so that the unknown parameter a is li nk ed 
to a random variable whose distribution is free of unknown parameters. We can achieve this by 

multiplying the left addend by 1, in the form of 


l (n-l)S 2 _ 

V(n~l)S 2 ' 


e 


Xn -X 


■ + Y-- 


-X ( n-l)s : 


{n-\)S : 


+ Y- 


-X (n-l)S 2 




+ Y. 


(5) 


(n -1)^ 2 

Now <J is linked to ---- which was noted in step 1 as having a parameter-free distribution. 

o~ 

S 2 

After V = (n - 1)—— is substituted in equation 5, all unknown parameters are removed from 0 : 


0 = 


Xn 


X \{n-l)S : 


Syfn^ I 


+ Y = 


Xp-X 

Syfn^l 


yfV +Y. 


( 6 ) 


Step 3 : At this point, we replace the remaining sufficient statistics in equation 6 with their 


observed values; that is, 


x 0 - x 
s-Jn -1 


yfV+Y. 


Step 4 : Finally, we rewrite the random variables V and Y in their original form (as functions of 
the sufficient statistics) to obtain the generalized pivotal quantity: 


r(x, S 2 ; x, s 2 ,ju,a) = 


x 0 - x j(n — l)iS 2 X — n 


i-Jn — 1 


( 7 ) 


As was demonstrated in the development of R , its distribution is free of any unknown parameters. 
Thus, condition 1 is satisfied. Furthermore, the observed value of R is 


r[x, s , x, s ,//,rr) 


X 0~ X S 2 ■ X ~M _n 

•\ 2 + 

.S' v a CJ 


thus satisfying condition 2. Therefore, R is a generalized pivotal quantity whose percentiles can 
be used to obtain a confidence interval for 6. Using the normal cumulative distribution, <!>(•), 

we then ultimately obtain a desired confidence interval for the probability of non-perforation. 


5 



3. Lower Confidence Bound (LCB) 


Since effective armor offers a high probability of non-perforation, engineers and management will 
most often be interested in a lower confidence bound, since this would represent a “pessimistic” 


bound on the true probability of non-perforation, <f> 


x 0 ~M 


A (l-a)l00% LCB for 6 = 


% 


is that value B, for which 


1 - a = P{B l < R) 


is satisfied. 4 We can estimate a generalized LCB for 0 by randomly generating a large number 
of observations of R and selecting a value for which equation 5 is empirically satisfied. For 
instance, with a = .05 , one could start by generating 100,000 observations of R, denoted by 
/ 1 , / 2 ,..., / 1 000 oo • If the ordered observations are written as /j l ( < n 2 ) < ... < n 100000 ), then a 95% 

LCB for 6 would be any real value B L so that /] 5000 ) < B L < >] 5m , j, such as B L = ^ 5000 ^ — ^ 5001 ^ . 

However, equation 8 can be expanded to get 


1 - a = PB L ± X -^S + ^ 

s \ a 2 (7 


P B l < 


x n -x V 



+ Y\. 


Notice that the random variable Y in equation 9 is equal in probability to W r , in which Z is a 

/ v« 


standard nonnal random variable (having mean 0 and variance 1). By using this fact and 
rearranging the terms in the last probability statement, we have 


i- a= J*!zl]JL + z kBl 

[ s \n -1 


x, -x 



sfv + Z ^ yfiiB , 


Z-^fnB, ^ x — x () 

= P — . L > -- 

s 


pT 

n—l,—'JnB L a / 


4 It is implicit that the upper bound equals 1 in the probability statement of equation 8, i.e., 1 - a = P(B L <R< 1). 
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In equation 10, T is a non-central t random variable with n -1 degrees of freedom and 


non-centrality parameter - 4nB L (see for example, Casella and Berger, 1990). However, a non¬ 
central t random variable with non-centrality parameter - 4nB L is the mirror image of a non¬ 
central t random variable with non-centrality parameter 4nB L (Johnson and Kotz, 1970). 
Therefore, 


\-a = 


P 



( 11 ) 


The final probability expression of equation 11 is the CDF of a non-central t random variable 
with n — 1 degrees of freedom and non-centrality parameter 4nB L , i.e., 

f \ 


1 -a = 




( 12 ) 


A lower confidence bound for G = 


_x 0 -ju 


is B l , the solution to equation 12. When we exploit the 


strict monotonicity of the normal distribution function, a (l -«)l00% LCB for the probability of 


non-perforation, <T> 


v O ) 


IS 


*A). 


4. Upper Confidence Bound (UCB) 


A UCB would represent an “optimistic” bound on the true probability of non-perforation. As such, 
it is of little practical value and rarely calculated. However, its derivation is briefly discussed here 
for completeness and as a precursor to the development of a two-sided confidence interval. 


Following a similar progression to that of the previous section, a (l -«)l00% UCB for 


x — (1 

0 - ——— is that value B r , satisfying 


f 


1 -a = P(R<B L ,) = P 




T ^ 

n-\,pnB v ^ / 


V 




Therefore, 


a = G 


n-\,PnB, 


f \ 

X 0 -X 




(13) 


V / vn J 
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( 

and a (l-«)lOO% UCB for the probability of non-perforation, <£> 

v 


O ) 


is <&(£„). 


5. Two-Sided Confidence Interval 


A two-sided confidence interval might be of interest to the researcher simultaneously wanting 
lower and upper limits for the probability of non-perforation, neither of which is automatically 
set to their extreme value of 0 or 1, respectively. 

A (l-«)l00% two-sided confidence interval for d = ~——— is given by values C L and C v 

<J 

satisfying 

^ = P(R<C L )mi^ = P(R>C v ) 


Following steps similar to those used in determining one-sided confidence bounds, two-sided 


confidence limits for 6 - 


_Xq~M 


are solutions to 


l-- = G r c 

2 n-\,yjnC L 


r 


x n - X 


V 




(14) 


and 


f \ 
x 0 -x 


— = G 

2 n—1,-JnCfj 

V./ J 

Thus, a (l -a)l00% two-sided confidence interval for the probability of non-perforation, 


(15) 


<f> 


( x 0 -/C 

l <7 , 


,is ^(cj,^)). 


6. Iterative Solutions for the Confidence Bounds and Confidence Limits 


As an example, consider solving equation 12 for the LCB of 0 = 


a 


the equation so that the right side equals 0: 


We start by rewriting 
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In equation 16, n, x 0 , x , s, and a are fixed constants, and the only unknown is B L . Therefore, 


if we let H(B L ;n,x 0 ,x,s,a) = 


+ a- 1, the problem becomes one of solving for 


the root of the function H. We can show if to be a monotonic decreasing function in B L by 
recognizing that the non-central t distribution function is monotonic decreasing in the non¬ 


centrality parameter (see figure 1) or by formal proof (see appendix A). Therefore, one can use 
the bisection method or other root-finding algorithm to solve for B L . Finally, the LCB for the 
probability of non-perforation is ( P{B i ). 


The bisection method can be easily programmed into most mathematical software packages. 
MATLAB 5 and Mathematica 6 programs appear in appendices B through G for calculating an 

LCB, a UCB, and a two-sided confidence interval for 




Figure 1. Two non-central t density functions, f and f 2 , with 


respective non-centrality parameters \ and /l,. (As the non¬ 
centrality parameter increases, the area to the right of X 0 
under the non-central t curve decreases.) 


5 MATLAB is a registered trademark of The MathWorks. 

6 Mathematica is a registered trademark of Wolfram Research, Inc. 
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7. Confidence Bands 


By plotting the function <J> 




for all real numbers x, one obtains a smooth, monotonically 


v * ) 

increasing function that is a simple estimate for the entire CDF. A lower confidence band for the 
CDF is a plot of the function <P(B l (x)) . This band lies below the estimated CDF and enables one 


to see the region above the band in which the entire actual CDF lies (with the stated level of 
confidence). If one fixes the probability level on the vertical axis, then the corresponding x-axis 
value from the lower confidence band is an upper tolerance bound. In the context of armor 
design, the upper tolerance bound is very important as it states with a specified level of 
confidence what annor thickness is needed to stop a (high) percentage of projectiles from 
perforating the materiel. A MATLAB program for generating a lower confidence band appears in 
appendix H. 


Upper confidence bands and two-sided confidence bands are constructed in a similar fashion; 
however, they are likely to be of much less interest to the annor designer. Although MATLAB 
programs for their construction are not included here, the interested reader could easily tailor the 
program in appendix H to accomplish this. 


8. An Application 


Table 1 gives the penetrations of 14 projectiles into an extended armor pack. Engineers plan to use 
115 units as the armor thickness. To determine if this thickness will provide enough protection, an 
estimate is desired for the probability that the next projectile fired will not penetrate deeper than 
115 units. Letting X be the penetration of this next projectile, we seek an estimate for 
p(x < 115). 


Table 1. Penetration depths into armor (no units specified). 


47 

81 

89 

99 

114 

59 

86 

90 

100 

118 

80 

88 

90 

113 



The data are first checked for normality by Lilliefors Test (see Conover, 1980). The value of the 
test statistic is 0.1713, which corresponds to a P-value greater than 0.20. Therefore, the 
assumption of normality is not rejected. 

Summary statistics for the sample data are x = 89.571 and .v 2 = 390.725. Equation 2 is used to 
obtain a point estimate for P(X < 115). The solution is 


10 




115-89.571 


<*>(1.286)= 0.901. 


P(X < 115) = 4> 


V390.725 


To achieve a 95% LCB for P(X <115) , we start by using equation 12: 

7 13,VI 4B l 


.95 = G„ rrTr. (4.813). 


Using the bisection method, we obtain a solution of B L = 0.6672. Therefore, a 95% LCB for 
P(X <115) is <f>(0.6672) = .7477. That is, for an annor of thickness 115 units, one can be 95% 
confident that the probability of non-perforation is at least 74.77%. 

At this point, a reasonable question might be “With 95% confidence, what armor thickness 
would offer at least 90% protection against perforation?” The lower confidence band for ), 

shown in figure 2, will help determine this. We begin by extending a horizontal (dashed) line 
from the “y”-axis at P(non - perforation) = .90 to the lower confidence band and then dropping 

to the x-axis to find the desired armor thickness. This thickness of 131.25 is the 90% upper 
tolerance bound. 


1 

0.9 

0.8 

z 0.7 
O 

0.6 
0.5 

LLI 
Q. 

± 0.4 
O 
z 

0- 0.3 
0.2 
0.1 
0 

20 40 60 80 100 120 140 160 180 

ARMOR THICKNESS 

Figure 2. Estimated probability of non-perforation and 95% LCB as a function of armor 
thickness. (The line formed by the bounds is referred to as the lower confidence 
band. From the band, one is able to determine what armor thickness yields a certain 
minimum probability of non-perforation, e.g., with 95% confidence, one can state 
that a thickness of 131.25 will stop at least 90% of projectiles from perforating the 
material.) 


h- 

£ 

o 

LL_ 

cc 
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A 95% two-sided confidence interval for P(x <115) is calculated from equations 12 and 13. It 


oc 

has a lower limit given by $>{C L ) in which C, is the solution to 1 - — = G ;1 j nC 

V / Vn J 

.975 = G , (4.813). The upper limit is given by f h(C ( ) in which C u is the solution to 


x 0 -x 


or 


13,Vl4C £ 
f 

- = G 

2 n-l,PnC u 


Xn ~ X 


or .025 = G 13 (4.813). Using the bisection method again, we obtain 


V / J 

C L = .5566 and C v =1.9901. Therefore, the lower confidence limit is d>(.5566) = .7111, and the 
upper confidence limit is <£>(l.990l) = .9767 . We can be 95% confident that the true probability 
of non-perforation is between 71.11% and 97.67%. 


9. Summary 


Reporting an interval estimate along with a point estimate gives designers a sense of the error of 
estimation which is a function of the sample size used in the study and the inherent variability of 
each projectile-armor interaction. Using a generalized pivotal quantity, we have achieved an 
interval estimate for the probability that a homogeneous armor plate of specified thickness will 
successfully stop a projectile from perforating. In the derivation of this interval estimate, we 
have assumed that the penetration depths are random and that they follow a normal distribution 
of unknown mean and variance. 

The interval limit (or limits if a two-sided interval is desired) is a function of the sample mean 
and sample standard deviation of penetration depths, the number of shots fired, the thickness of 
the armor and percentiles of a non-central t distribution. However, since the non-centrality 
parameter associated with this distribution is a function of the interval limit, numerical methods 
are required to obtain the final solution. MATLAB and Mathematica codes to perform these 
calculations are provided in appendices B through G. 

Confidence bands that graphically display the entire the relationship between armor width and 
bounds on the probability of non-perforation follow naturally. These bands may be of help to 
armor designers in selecting an armor width that will provide a high degree of protection against 
specific enemy threats. 
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Appendix A. Proof of the Monotonicity of the Function H 


Show that H( ) is a monotonic decreasing function. That is, show that if A l < Z ,, then 

h{a 1 )>h{a 2 ). 

Proof: 

Let n, x 0 , x ,s, a, A l , and A, be fixed constants with A l < Z, . Let X x and The independent 
random variables, X t being normally distributed with mean A l and variance 1, and Y following a 
chi-squared distribution with n -1 degrees of freedom. The support of X l is the set of all real 
numbers, and the support of Y is the set of all positive reals. 


Now define W ■ 


yiY x — x 

-——— . The support of W is either the set of all positive real numbers (if 

n -1 .s' 


x 0 > x ) or the set of all non-positive real numbers (if x 0 <x). However, in general, we denote 
the support of W by 0 . 

P(X j < W) can be expressed as a double integral of the joint density of X l and W, namely, 


w 

p(x x <w)= IJa , (a')/ h , (w jt/x dw. 


(A-l) 


The double integral of equation A-l, can be rewritten as the sum of two double integrals: 

W-(A2 ~^\ ) w 

p{x, < W ) = J j fx x U)f w ( w)dx t/vv + J J f Xi (x)f v ( w)dx dw. 

0 —OO 0 W—(Aq —Ay ) 


(A-2) 


Each of these double integrals is positive, so by dropping the latter of the two, we obtain the 
inequality, 


W -(^2 ~ k \ ) 

P[X j < W) > J \f Xi {x)f w (ujdx dw . 


(A-3) 


However, the double integral in equation A-3, represents a probability, namely, 

P(X j < W - (A 2 - A x )) , which in turn equals P{X l + A 2 - A 1 < W) . So, 

P{X l <W)> P(X x +A 2 -/L, < W). (A-4) 

Now define X 2 = X x + (A 2 - A,) ; then X 2 is normally distributed with mean A 2 and variance 1. 
Therefore, 


p(x 1 <w)>p(x 2 <w). 


(A-5) 
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Consider the standard normal random variable, Z, having mean 0 and variance 1. Employing the 
equalities X, = Z + A, and X 2 = Z + A-, and the definition of W, we expand the above result to 

obtain 


P\ Z + /L < 


nY x ( . 



>P\ Z + /L < 


nY x n 



(A-6) 


Rearranging the expressions inside each of the probability statements, 



By the definition of a non-central t random variable, 



(A-7) 


(A-8) 


Each of these probabilities is the cumulative non-central t distribution function evaluated at 




Next, we add the constant a -1 to both sides of the inequality: 



Using the definition of the function //(■), we obtain the desired result. 

h{\)> H{X 2 ) . 


(A-10) 


(A-11) 
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Appendix B. MATLAB Code for Probability of Non-perforation LCB 


MATLAB Release 14, Version 7, code for calculating an LCB for the probability of non¬ 




perforation, $1 ' 0 _ / “ |. The user may supply his or her own data and values for alpha, xO, and 
epsilon. 


% Define variables % 

alpha=.05; x0=115; epsilon=.000001; 

% Declare data % 

data=[47 59 80 81 86 88 89 90 90 99 100 113 114 118]'; 


% Calculate summary statistics and define K % 
n=length(data); xbar=mean(data); sigma=std(data); 

K=sqrt(n)*(x0- xbar)/sigma; 

% Initialize flag and bounds on non-centrality parameter % 
flag=0; lononcent=xbar-10*sigma; hinoncent=xbar+10*sigma; 

% Execute bisection method % 
noncent=(lononcent+hinoncent)/2; 
while flag==0 

if hinoncent-lononcent < epsilon 
flag=l; 

elseif nctcdf(K,n-1,noncent)+alpha-l < 0 

hinoncent=noncent; noncent=(noncent+lononcent)/2; 

else 

lononcent=noncent; noncent=(noncent+hinoncent)/2; 

end 

end 


%Print lower confidence bound for probability of non-perforation % 
fprintf('%3.Of%% LCB for P(non-perforation) = %6.4f\n',... 

100* (1-alpha),normcdf(noncent/sqrt(n)) ) 
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Intentionally left blank 
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Appendix C. Mathematica Code for Probability of Non-perforation LCB 


Mathematica Version 5.1 code for calculating an LCB for the probability of non-perforation. 




X — // 

—- . The user may supply his or her own data and values for alpha, xO, and epsilon. 

V O' ) 


(* Remove warning messages that do not affect results but clutter output *) 
Off[General: :spell]; Off[General: :spell1] ; 

Off[NIntegrate::slwcon]; Off[NIntegrate::ncvb]; 


(* Enable use of statistical tools *) 
<<Statisties'ContinuousDistributions' 
<<Statisties'NormaIDistribution' 


(* Clear variables for use *) 

Clear[data, n, noncent, xO,a,epsilon,lononcent,hinoncent,flag,K,ndist] 

(* Define variables *) 

alpha=.05; x0=115; epsilon =.000001; 

(* Declare data *) 

data={47,59,80,81,86,88,89,90,90,99,100,113,114,118}; 

(* Calculate summary statistics and define K *) 

n=Length[data]; xbar=Mean[data]//N; sigma=StandardDeviation[data]//N; 
K=(xO-xbar)/(sigma/Sqrt[n])//N; 

(* Initialize flag and bounds on non-centrality parameter *) 
flag=0; lononcent=xbar-10sigma//N; hinoncent=xbar+10sigma//N; 

(* Execute bisection method *) 

H:=CDF[NoncentralStudentTDistribution[n-1,noncent],K]+alpha-l 
While[flag==0,noncent=Mean[{lononcent,hinoncent}]; 

If[hinoncent-lononcent<epsilon,flag=l. 

If[H<0, hinoncent=noncent,lononcent=noncent]]] 

(* Print lower confidence bound for probability of non-perforation *) 
Print[TraditionalForm[StringForm["''0% LCB for P(non-perforation) = 

100 (1-alpha),CDF[NormaIDistribution[0,1],noncent/Sqrt[n]]]]] 


19 




Intentionally left blank 
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Appendix D. MATLAB Code for Probability of Non-perforation UCB 


MATLAB Release 14, Version 7, code for calculating a UCB for the probability of non- 

( jp — n b 

perforation, <f> —- . The user may supply his or her own data and values for alpha, xO, and 

\ a ) 

epsilon. 

% Define variables % 

alpha=.05; x0=115; epsilon=.000001; 

% Declare data % 

data=[47 59 80 81 86 88 89 90 90 99 100 113 114 118]'; 


% Calculate summary statistics and define K % 
n=length(data); xbar=mean(data); sigma=std(data); 
K=sqrt(n)*(x0- xbar)/sigma; 


% Initialize flag and bounds on non-centrality parameter % 
flag=0; lononcent=xbar-10*sigma; hinoncent=xbar+10*sigma; 


% Execute bisection method % 
noncent=(lononcent+hinoncent) / 2 ; 
while flag==0 

if hinoncent-lononcent < epsilon 
flag=l; 

elseif nctcdf(K,n-1,noncent)-alpha < 0 

hinoncent=noncent; noncent=(noncent+lononcent)/2; 

else 

lononcent=noncent; noncent=(noncent+hinoncent)/2; 

end 

end 


% Print upper confidence bound for probability of non-perforation % 
fprintf('%3.Of%% UCB for P(non-perforation) = %6.4f\n',... 

100* (1-alpha),normcdf(noncent/sqrt(n))) 


21 




Intentionally left blank 


22 



Appendix E. Mathematica Code for Probability of Non-perforation UCB 


Mathematica Version 5.1 code for calculating a UCB for the probability of non-perforation. 




X — // 

—- . The user may supply his or her own data and values for alpha, xO, and epsilon. 

V O' ) 


(* Remove warning messages that do not affect results but clutter output *) 
Off[General: :spell]; Off[General: :spell1] ; 

Off[NIntegrate::slwcon]; Off[NIntegrate::ncvb]; 


(* Enable use of statistical tools *) 
<<Statisties'ContinuousDistributions' 
<<Statisties'NormaIDistribution' 


(* Clear variables for use *) 

Clear[data, n, noncent, xO,a,epsilon,lononcent,hinoncent,flag,K,ndist] 

(* Define variables *) 

alpha=.05; x0=115; epsilon =.000001; 

(* Declare data *) 

data={47,59,80,81,86,88,89,90,90,99,100,113,114,118}; 

(* Calculate summary statistics and define K *) 

n=Length[data]; xbar=Mean[data]//N; sigma=StandardDeviation[data]//N; 
K=(xO-xbar)/(sigma/Sqrt[n])//N; 

(* Initialize flag and bounds on non-centrality parameter *) 
flag=0; lononcent=xbar-10sigma//N; hinoncent=xbar+10sigma//N; 

(* Execute bisection method *) 

H:=CDF[NoncentralStudentTDistribution[n-1,noncent],K]-alpha 
While[flag==0,noncent=Mean[{lononcent,hinoncent}]; 

If[hinoncent-lononcent<epsilon,flag=l. 

If[H<0, hinoncent=noncent,lononcent=noncent]]] 

(* Print upper confidence bound for probability of non-perforation *) 
Print[TraditionalForm[StringForm["''0% UCB for P(non-perforation) = 

100 (1-alpha),CDF[NormaIDistribution[0,1],noncent/Sqrt[n]]]]] 
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Intentionally left blank 
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Appendix F. MATLAB Code for Probability of Perforation Two-Sided 
Confidence Interval 

MATLAB Release 14, Version 7, code for calculating a two-sided confidence interval for the 

f % — 

probability of non-perforation, <£> ——— . The user may supply his or her own data and values 

for alpha, xO, and epsilon. 

% Define variables % 

alpha=.05; x0=115; epsilon=.000001; 

% Declare data % 

data=[47 59 80 81 86 88 89 90 90 99 100 113 114 118]'; 

% Calculate summary statistics and define K % 
n=length(data); xbar=mean(data); sigma=std(data); 

K=sqrt(n)*(x0- xbar)/sigma; 

% Initialize flag and bounds on non-centrality parameter % 
flag=0; lononcent=xbar-10*sigma; hinoncent= xbar+10*sigma; 

% Execute bisection method for lower limit % 
noncent=(lononcent+hinoncent)/2; 
while flag==0 

if hinoncent-lononcent < epsilon 
flag=l; 

elseif nctcdf(K,n-1,noncent)+alpha/2-l < 0 

hinoncent=noncent; noncent=(noncent+lononcent)/2; 

else 

lononcent=noncent; noncent=(noncent+hinoncent)/2; 

end 

end 

% Print lower confidence limit for probability of non-perforation % 
fprintf('%4.Of%% Cl for P(non-perforation) = (%6.4f, 

100*(1-alpha) , normcdf(noncent/sqrt(n))) 

% Re-initialize flag and bounds on non-centrality parameter % 
flag=0; lononcent=xbar-10*sigma; hinoncent= xbar+10*sigma; 
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% Execute bisection method for upper limit % 
noncent=(lononcent+hinoncent)/2; 
while flag==0 

if hinoncent-lononcent < epsilon 
flag=l; 

elseif nctcdf(K,n-1,noncent)-alpha/2 < 0 

hinoncent=noncent; noncent=(noncent+lononcent)/2; 

else 

lononcent=noncent; noncent=(noncent+hinoncent)/2; 

end 

end 

% Print upper confidence limit for probability of non-perforation 
fprintf('%6.4f)\n',normcdf(noncent/sqrt(n))) 
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Appendix G. Mathematica Code for Probability of Perforation Two-Sided 
Confidence Interval 


Mathematica Version 5.1 code for calculating a two-sided confidence interval for the probability 


of non-perforation, <f> 
xO, and epsilon. 


f *o-/^ 


. The user may supply his or her own data and values for alpha, 


(* Remove warning messages that do not affect results but clutter output *) 
Off[General::spell]; Off[General::spell1]; 

Off[NIntegrate::slwcon]; Off[NIntegrate::ncvb]; 


(* Enable use of statistical tools *) 
<<Statisties'ContinuousDistributions' 
<<Statisties'NormaIDistribution' 


(* Clear variables for use *) 

Clear[data,n,noncent,xO,a,epsilon,lononcent,hinoncent, flag, K, 
ndist,lowerlimit,upperlimit] 

(* Define variables *) 

alpha=.05; x0=115; epsilon =.000001; 

(* Declare data *) 

data={47,59, 80,81,86,88,89,90,90,99,100,113,114,118}; 

(* Calculate summary statistics and define K *) 

n=Length[data]; xbar=Mean[data]//N; sigma=StandardDeviation[data]//N; 
K=(xO-xbar)/(sigma/Sqrt[n])//N; 

(* Initialize flag and bounds on non-centrality parameter *) 
flag=0; lononcent=xbar-10sigma//N; hinoncent=xbar+10sigma//N; 

(* Execute bisection method for lower limit *) 

H:=CDF[NoncentralStudentTDistribution[n-1,noncent],K]+alpha/2-1 
While[flag==0, noncent=Mean[{lononcent,hinoncent}]; 

If[hinoncent-lononcent<epsilon,flag=l. 


27 




If[H<0,hinoncent=noncent,lononcent=noncent]]] 
lowerlimit = N[CDF[NormalDistribution[0,1],noncent/Sqrt[n]]]; 

(* Re-initialize flag and bounds on non-centrality parameter *) 
flag=0; lononcent=xbar-10sigma//N; hinoncent=xbar+10sigma//N; 

(* Execute bisection method for upper limit *) 

H:=CDF[NoncentralStudentTDistribution[n-1,noncent],K]-alpha/2 
While[flag==0,noncent=Mean[{lononcent,hinoncent}]; 

If[hinoncent-lononcent<epsilon,flag=l. 

If[H<0, hinoncent=noncent,lononcent=noncent]]] 
upperlimit = N[CDF[NormalDistribution[0,1],noncent/Sqrt[n]]]; 

(* Print confidence interval for probability of non-perforation *) 
Print[TraditionalForm[StringForm["''0% Cl for P(non-perforation) = 
100 (1-alpha),lowerlimit,upperlimit]]] 
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Appendix H. MATLAB Code for Probability of Non-perforation CDF and 
Lower Confidence Band 


MATLAB Release 14, Version 7, code for plotting an estimated CDF and lower confidence 
band. The user may supply his or her own data and values for alpha, xO, and epsilon. 

% Define variables % 
alpha=.05; epsilon=.0001; 

% Declare data % 

data=[47 59 80 81 86 88 89 90 90 99 100 113 114 118]'; 

% Calculate summary statistics % 

n=length(data); Xbar=mean(data); SDev=std(data); 

% Initialize counter for points along the X-axis % 
i=0; 

% Declare armor thickness % 
for xO=Xbar-3*SDev:SDev/10:Xbar+4*SDev 
i=i+l; x(i)=x0; 

% Calculate plug-in estimate for probability of non-perforation % 
pointest(i)=normcdf((xO-Xbar)/SDev) ; 

% Define K % 

K=sqrt(n)*(xO-Xbar)/SDev; 

% Initialize flag and bounds on non-centrality parameter % 
flag=0; lononcent=-1000; hinoncent=1000; 

% Execute bisection method for lower confidence bound % 
noncent= (lononcent+hinoncent)/2; 
while flag==0 

if hinoncent-lononcent < epsilon 
flag=l; 

elseif nctcdf(K,n-1,noncent)-1+alpha < 0 
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hinoncent=noncent; noncent=(noncent+lononcent)/2; 

else 

lononcent=noncent; noncent=(noncent+hinoncent)/2; 

end 

end 

lcb(i)=normcdf(noncent/sqrt(n)); 

end 

% as a function of armor thickness, plot non-perforation probability curve 
% in blue, and lower confidence band for probability of non-perforation in 
% red % 

plot(x,pointest, 'b-',x,lcb, 'r-' ) 

xlabel('ARMOR THICKNESS'); ylabel('P(NON-PERFORATION)') 

text(Xbar-.25*SDev,.5,'Estimated','Color','b','Rotation',69,... 

'HorizontalAlignment','Center') 

text(Xbar+.7 5*SDev, .5, 'LC Band', 'Color', 'r ' , 'Rotation',67, . . . 

'HorizontalAlignment','Center') 
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List of Symbols Used 


Scalar Values 

n number of observations in sample data 

r. the i' h observed value of the random variable R 

r (; , the i' h ordered, observed value of the random variable R 

s observed sample standard deviation 

x, the i' h observed value from a data set; in this paper, the observed depth of penetration 
from the i h projectile 
x 0 a specified thickness of armor 

x observed sample mean 

z p the value from the standard normal distribution having an area of p to its right under the 
density curve 

a level of significance 

o population standard deviation for a normally distributed random variable 

6 parameter for which a confidence interval is constructed; may actually be a function of 

one or more parameters 

j.i population mean for a normally distributed random variable 
£ a set of one or more nuisance parameters 

Random Variables 

R( ) generalized test variable; a function of random data, observed data, parameter of interest 
and perhaps other nuisance parameters 
S sample standard deviation of depths of penetration 

T student’s t random variable 

V chi-square random variable 

X depth of penetration for a randomly selected projectile 

X sample mean depth of penetration 
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Y a normal random variable with mean zero and variance — 

n 

Z a standard normal random variable, with mean zero and variance one 

Functions 

F{ ) cumulative distribution function 
F( ) estimated cumulative distribution function 
G( ) student’s t cumulative distribution function 
<f>( ) standard normal cumulative distribution function 

Miscellaneous 

B l lower confidence bound 

B li upper confidence bound 

C L lower confidence limit 

C lf upper confidence limit 

P( ) probability of the parenthesized expression 

P{ ) estimated probability of the parenthesized expression 
E summation 
A estimate of 
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